Intrinsic spin-orbit coupling in superconducting 5-doped SrTi03 heterostructures 
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We report the violation of the Pauli limit due to intrinsic spin-orbit coupling in SrTi03 het- 
erostructures. Via selective doping down to a few nanometers, a two-dimensional superconductor 
is formed, geometrically suppressing orbital pair-breaking. The spin-orbit scattering is exposed by 
the robust in-plane superconducting upper critical field, exceeding the Pauli limit by a factor of 
4. Transport scattering times several orders of magnitude higher than for conventional thin film 
superconductors enables a new regime to be entered, where spin-orbit coupling effects arise non- 
pert urbatively. 



Unconventional superconductivity is a subject of great 
theoretical and experimental interest A central is- 

sue in this field is the discovery and understanding of 
non-trivial pairing mechanisms, such as the spin-triplet 
Cooper pair, which has been explicitly investigated in 
heavy fermions [1 , Sr2Ru04 [4 , and crystals with broken 
inversion symmetry [5]. Recently, novel pairing has also 
been predicted in two-dimensional systems breaking in- 
version symmetry j6l [7] . Experimentally, measurements 
of the superconducting upper critical held H c2 give vital 
information. In particular, violations of the Pauli para- 
magnetic limit 8 , 9 can be used to unravel the nature of 
the electron spins in the superconducting state. Notably, 
the presence of spin-orbit coupling (SOC) can be quan- 
tified [10J . as demonstrated by the studies of metal 
thin-film superconductors and bilayer systems where 
interface SOC drastically enhances H C 2 [12] . 

Electron-doped SrTi0 3 (STO) has attracted much at- 
tention as the lowest-density superconductor [T3] with 
high-mobility [T3]. These characteristics enable the cre- 
ation of novel low dimensional systems |15j . and are 
vital to shed light on the rich physics present at the 
LaA10 3 /SrTi0 3 (LAO/STO) interface, where the pres- 
ence of the Rashba spin-orbit interaction has been dis- 
cussed, affecting both the normal and superconducting 
state transport properties [TBI HI] • However, despite the 
fact that the conduction band structure of STO is sim- 
ilar to p-type GaAs [HI [19] , the latter a model system 
for spintronics, the role of possible intrinsic SOC in the 
transport properties of doped STO is still unclear. 

In this Letter, we study the violation of the supercon- 
ducting Pauli limit due to intrinsic SOC in a systematic 
series of symmetric, doped STO heterostructures. Using 
the ^-doping technique, we selectively add Nb dopants in 
a narrow region inside an otherwise continuous undoped 
STO host crystal. As the thickness of the dopant layer is 
reduced, the destruction of superconductivity by orbital 
pair-breaking is geometrically suppressed, and the super- 
conducting H C 2 is enhanced for magnetic fields applied 
parallel to the dopant plane. In the thin regime, when 



the dopant layer is just a few nanometers thick, the su- 
perconductivity is robust beyond the conventional Pauli 
limit, demonstrating the presence of spin-orbit scatter- 
ing (SOS) in the STO. Moreover, due to the absence of 
a surface or interface close to the dopant plane, and the 
spreading of the electron wavefunctions into the undoped 
STO, the electronic mean free path does not collapse as 
the dopant thickness decreases. Thus we preserve trans- 
port scattering times several orders of magnitude higher 
than for conventional thin film superconductors. In this 
regime, the intrinsic band SOC effects arise as a non- 
perturbative correction to the transport, despite the rel- 
atively long absolute SOS times. 

The samples were fabricated with various thicknesses 
of 1 at. % doped Nb:SrTi0 3 (NSTO) films embedded 
between cap and buffer layers of undoped STO, using 
pulsed layer deposition. High-temperature growth, above 
~ 1050 °C, in a low oxygen partial pressure of less than 
10~ 7 Torr was chosen to achieve high-quality STO films, 
by managing the defect chemistry of the strontium and 
oxygen vacancies [3U]. On a Ti02 terminated STO (100) 
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FIG. 1. (color online) (a) Sheet resistance R, normalized by 
the normal-state value _R n as a function of temperature T. 
Numbers refer to the <5-doped layer thickness d in nm. (b) 
Superconducting transition temperature T c versus d. T c is 
defined by the temperature at the half value of the normal- 
state resistance; 10 %-90 % width of resistance is shown as 
an error bar. 
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FIG. 2. (color online) (a) Angular dependence of the upper critical field H C 2 at 50 mK, normalized by the value at 9 = 
°. Dashed curves are fitting results obtained to Tinkham's model. Results for representative samples (numbers refer to d in 
nm) and a bulk f at. % NSTO substrate (Subs.) are shown, (b) Normalized perpendicular upper critical field H^/ '//^(O 
K) plotted as a function of the reduced temperature t = T/T c . H^(0 K) was obtained by extrapolation to T = K from a 
fitting to the data in (b) over the range of 0.7 < t < 1. (c) Normalized parallel upper critical field with same fitting procedure, 
(d) Ginzburg-Landau coherence length £gl (open rectangles) and dTinkham (closed diamonds) versus d. Dashed line is d = 

C?Tinkham- 



substrate, a 100 nm undoped STO buffer layer was first 
grown, followed by the 1 at. % NSTO layer with various 
thicknesses in the range 3.9 nm < d < 457 nm. A 100 
nm undoped STO cap layer was grown above the doped 
layer, to prevent surface depletion |21j . Post-annealing in 
a moderate oxidizing condition was used to fill oxygen va- 
cancies formed during growth. Transport measurements 
were made using a standard four-probe method with sam- 
ple cooling achieved using a dilution refrigerator with an 
in-situ rotator. For zero field measurements, the residual 
magnetic field was reduced below an absolute value of 
(IqH = 0.1 mT, where [1q is the vacuum permeability. 

All samples were superconducting at low temperatures, 
as shown in Fig. [I] (a). The transition temperatures T c 
defined by the temperature below which the resistance 
was 50 % of the normal-state value, were in the range 
253 mK < T c < 374 mK, as shown in Fig. [I] (b). All 
samples, except for the two thinnest, showed sharp 10 %- 
90 % transition widths (~ 10 mK). While samples with 
thickness d > 8.8 nm showed relatively constant T c (~ 
260 mK), several thinner samples showed a higher T c 
while maintaining a sharp transition, suggesting possible 
changes to the superconducting properties close to the 
two-dimensional (2D) limit. We note that although the 
transition broadening in some of the thinner samples may 
relate to inhomogeneities, this is also reminiscent of the 
suggested Bose metal phase between the superconducting 
and insulating states [2"2"| . 

Firstly, the anisotropy of H C 2 was used to measure the 
dimensionality of the superconductivity. We investigated 
the variation of H C 2 by rotating the sample with respect 
to the magnetic field, as shown in Fig. [2] (a). A bulk 1 at. 
% NSTO substrate was also measured as a reference. As 
d decreased, a clear modulation of H c2 as a function of the 
angle 9 between the magnetic field and the sample plane 
was found. Here H C 2 was defined as the field at which the 
resistance was half that of the normal state. For samples 



with d < 99 nm, excellent fits to these data could be 
made using Tinkham's model |23j . which is valid when 
the superconducting thickness is less than the Ginzburg- 
Landau coherence length dxinkham < Cgl(0) [21] ■ These 
fits are shown in Fig. [2] (a) . 

The dimensional crossover of superconductivity is more 
clearly demonstrated by the temperature dependence of 
H c2 , therefore we next measured H^(t) and H^ 2 (t), the 
out-of-plane (8 = 90 °) and in-plane (9 = 0°) upper criti- 
cal fields respectively (here t = T/T c ), as shown in Figs. [2] 
(b) and (c) . In the perpendicular field geometry, all sam- 
ples showed a linear temperature dependence. In the 
parallel field geometry, for d < 99 nm, however, H\ 2 {t) 
showed a clear square root form, which is characteris- 
tic of the 2D superconducting state. These data clearly 
demonstrate a three-dimensional (3D) to 2D crossover 
of the superconducting character as a function of d. By 
estimating c?Tinkham and £gl(0) from the and H^ 2 
data, we find that G?Tinkham decreases in proportion to 
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FIG. 3. (color online) Hj: 2 (circles), (squares, scaled by 
factor of five), and the Pauli paramagnetic limit H£ (trian- 
gles) plotted by 1/d. H^ 2 and are 50 mK data. Dashed 
lines are guides to the eye. 
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FIG. 4. (color online) (a) Contour plot showing the deviation between the WHH simulation and the experiment by using H c i{0) 
data for d — 5.5 nm. The minimum point is located at A = 6.4, a — 0.14. (b) H C 2{6) of 5.5 nm thick sample at 50 mK. Dotted 
line is the best fit obtained from the WHH simulation. Inset: H^{T) and H^ 2 (T) data and the WHH theory fit (dotted line), 
(c) Variation of r so and r tr with d for the five thinnest samples, obtained from best fits to the WHH theory. An error bar of r so 
is given by assuming 10 % thickness variation of the superconducting layer. The ratio r tr /T so is also shown on the right axis. 



the growth thickness d, and in the thinnest sample is 
much smaller than £gl(0) ~ 100 nm, as plotted in Fig. [2] 
(d), confirming the 2D nature of the superconductivity. 

A crucial and intriguing aspect of the H} 2 (t) data 
is the violation of the Pauli paramagnetic limit. The 
Pauli paramagnetic limiting field [SI IS] is given by = 
Ao/\/2/iB, where /iq is the Bohr magneton (with a g- 
factor of two), and Ao = 1.76 fee ?c is the BCS supercon- 
ducting gap for a weak-coupling superconductor, where 
&b is Boltzmann's constant. This limit is appropriate, 
since via tunneling bulk doped STO is known to be in 
the weak-coupling regime . The variation of , H^ 2 , 
and if P as a function of 1/d is shown in Fig. [3J i?J 2 ex- 
ceeds the Pauli limiting field iJ P by a factor of more than 
four in the thinnest sample, while and Hf? remained 
essentially constant. 

In the case of a 2D superconductor in a parallel mag- 
netic field, if the sample is thin enough that orbital de- 
pairing is suppressed, spin paramagnetism is the domi- 
nant mechanism for destroying superconductivity |11) . In 
the presence of SOS, however, h\ 2 can be robust beyond 
the Pauli limit. It should be noted that the renormaliza- 
tion of normal-state properties by many-body effects |26j 
can also enhance the Pauli limit. However, Shubnikov-de 
Haas oscillations in these 5-doped samples showed that 
the electron mass is consistent with the band structure 
at low temperature [15J . Given the low electron-density, 
strong correlation effects near half-filling are also absent. 

To investigate the effects of SOS in more detail, we 
performed a numerical fit of the H c2 (6,t) data to the 
Werthamer-Helfand-Hohenberg (WHH) theory [TU] tak- 
ing into account corrections for the thin film case [24] , 
Within this theory, the two crucial fitting parameters, 
the orbital depairing parameter a and the SOS rate A so , 
are given by 
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37rfc B JcT S o 

where H is the Plank constant divided by 2ir, m* is the 
effective electron mass, D = VpT tr /3 is the diffusion con- 
stant, and vf is the Fermi velocity. r tr and t so are, re- 
spectively, the transport and SOS times. For various a 
and A so , we calculated the sum of the squares of the dif- 
ferences between the WHH model and the H c2 (9) data 
(we denote this sum as £), for the case d = 5.5 nm, as 
shown in Fig. [4] (a) . A unique minimum value of S was 
found, giving an excellent fit to the experimental data, 
as shown in Fig. @ (b). We obtain T tr = 6.2 x 10 s 
and t so = 6.8 x 10~ 13 s, from a and A so , respectively. 
This clearly indicates that the SOS is a highly significant 
factor affecting the normal and superconducting trans- 
port properties, in spite of its relatively small absolute 
value. This is due to an appealing point of the <5-doped 
STO structures, where Tt r is several orders of magnitudes 
higher than for conventional metal superconductors [57] . 

By performing similar analysis on data for the other 
samples where £/J 2 exceeded 7J P , we obtain the thickness 
dependences of t so and T tr , as shown in Fig.[I](c). A clear 
decrease of r so with decreasing d is found, while r tr is rel- 
atively unchanged. This dependence suggests that the 
SOS is not dominated by either the Elliott- Yafet mecha- 
nism where r so cx r tr , or the D'yakonov-Perel' mechanism 
where r so cx 1/ Tt r [28] . This rather unexpected result sug- 
gests that the SOS observed has a different origin. We 
next clarify this point by comparison with other systems. 

We emphasize that this combination of SOC with high 
mobility conduction electrons places our system is in a 
different regime compared to other thin film structures 
that violate the Pauli limit. For example, the use of 
heavy atoms to induce SOC in superconducting bilayers 
has been studied [12]. However, in this case, as is usual 
for conventional superconducting thin films, the mean 
free path collapses in the thin limit. A similar collapse 
occurs with substrate gating at the LAO/STO heteroin- 
terface |29j where an asymmetric confining potential and 
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Rashba SOC is expected [THl E]- <5-doping is a crucial 
determinant for this difference: since there is no obvious 
surface or interface surrounding the conducting layer, the 
scattering length due to disorder is unchanged (and even 
increased) with decreasing d [30] . Additionally, the sym- 
metry of the structure, giving rise to zero net effective 
electric field, means that Rashba SOC is absent. 

We can thus interpret this as intrinsic SOC of STO, 
due to the d-orbitals of the Ti atoms [18]. Indeed the 
bulk conduction bands of STO have a similar structure 
to the valence bands of GaAs, where non-perturbatively 
large SOC has been demonstrated [T0J- However in the 
case of STO there are few studies of electron SOS. It 
should be noted that calculations indicate that these con- 
fined (5-doped samples have a multiple subband structure 
[T5] . therefore the change of r so observed may relate to 
intersubband- induced spin-orbit interaction [31] . or in- 
tersubband scattering [32j . which in turn are influenced 
by changes of the band structure as d changes. 

The fact that the energy scale of the observed SOC (^2 
meV) is bigger than the superconducting gap (^40 ^eV) 
suggests the possibility of mixed spin-triplet and singlet 
states, giving rise to novel superconducting states in these 
low dimensional layers [51ITll33|. Moreover, in the normal 
state, the combination of high mobility conduction and 
SOC imply that STO can be usefully employed in a wide 
range of controllable spintronic architectures, including 
the spin Hall effect, which have until now been domi- 
nated by metals and traditional semiconductors. Thus 
this system is ideal for future measurements, which can 
be made over a range of densities both via growth con- 
trol [33], and field-effect gating, the latter simultaneously 
introducing Rashba contributions to the SOC [55] . 
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Tinkham's model 

The magnetic-field response of a two-dimensional su- 
perconductor can be described by Tinkham's model [T], 
which assumes that the thickness of the superconductor 
is thinner than the Ginzburg-Landau coherence length, 
^Tinkham < £gl- It should be noted that the model does 
not include the effects of spin-orbit scattering or the Pauli 
paramagnetic limit, and assumes an isotropic supercon- 
ducting wavefunction. According to the model, the an- 
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gular dependence of the upper critical field can be shown 
to be 

2 



H c2 (6) sm9 



H c2 (6)coat 



H. 



= 1, 
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c2 



where 8 is the angle between the magnetic field and the 
sample plane. The temperature dependence of the upper 
critical field in perpendicular and parallel field geometry 
is given by 



g«a(*) = 9 /", nV , (1-*). 
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where t = T/T c is the reduced temperature, <E>o = h/2e = 
2.07 x 1CT 15 Wb is the flux quantum, and £gl(0) is the 
Ginzburg-Landau coherence length extrapolated to T — 
K. From Eqs. [2] and [3j dTinkham and £gl(0) can be 
found using 
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(5) 



Thus, ^Tinkham can be calculated by measurement of 
and i?| 2 of a sample experimentally. As noted by Ben 
Shalom et al. [2], in the case of the LaAlOa/SrTiOa in- 
terface, the value of ^Tinkham is an upper bound on the 
thickness. In our case we find good agreement between 
the grown dopant layer thickness d and dxinkham in thick 
samples, but ^Tinkham deviates slightly from d as ijj 2 
exceeds the Pauli limiting field in the thinnest samples 
(d < 8.8 nm), indicating a limit of the model. 
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where D is the diffusion constant, r so the spin-orbit scat- 
tering time, and m the electron mass. To include the 
effect of the finite thickness of a thin film, the term h in 
Eq. [6] should be replaced by h ans (8), which is given by 
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where d is thickness of superconducting layer. 

In the case of a — A so = (no spin-orbit coupling, 
but also no Pauli paramagnetic limit), the above formula 
is reduced to the orbital term only, where H c2 is the 
solution of 

lnt + V,Q + ^Q)=0. (11) 

In fitting the data, the orbital only case (Eq. [TTj ) could 
not accurately fit the data for samples with H^ 2 larger 
than the Pauli paramagnetic limit. However we could 
obtain a very successful fit with the full WHH theory 
(Eq. [6]), as discussed and shown in the main text (Fig. 
4). 

Several sources of error should be considered in this 
fit. Firstly, it has been argued that cooling supercon- 
ducting ultra-thin films below 60 mK is extremely dif- 
ficult [B], thus the H c2 (t) data may artificially saturate 
at low temperatures due to a lack of cooling. However, 



WHH theory 

The Werthamer-Helfand-Hohenberg (WHH) theory [3] 
was used to more quantitatively fit the superconducting 
upper critical field, H c2 , data in the main text, in order 
to determine the spin-orbit scattering time in the system. 
Within this theory H c2 is the implicit solution of the 
equation 
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where ip is the digamma function. With a slight cor- 
rection |4 from the original WHH paper, the terms are 
defined as 

r DeH c2 



nk B T c 
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FIG. SI. Optimal value of A so from the WHH fitting for the 
d = 5.5 nm sample, depending on the thickness d used in the 
model. 
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FIG. S2. Variation of Hj} 2 with \ ao using a non-perturbative 
spin-orbit coupling model, compared to the original WHH 
model, adapted from Schopohl et al. [8|. 



by fitting to data at temperatures only above 75 mK, we 
find an increase of the value of t so of only 2 % compared 
to the full fitting curve. Secondly, the influence of error 
in the value of d used in Eq. [10] can also be considered. 



This effect is demonstrated in Fig. SI where the value 
of t so obtained from the fitting is plotted against d for 
the d = 5.5 nm sample. As is clear, sensitivity of the fit 
to variation in d gives rise to variation in r so ; we obtain 
t so = (6.8 ± 1.0) x KT 13 s for d = 5.5 ± 0.5 nm. 

We note that the WHH theory makes various simpli- 
fying assumptions: the superconductor should be in the 
dirty limit, where the electron mean free path is shorter 
than the BCS coherence length £ <C £bcs- Secondly, the 
spin-orbit scattering time is less than the total scatter- 
ing time t so <C 7"tr- In order to estimate these parame- 
ters, we assume a single band approximation with spher- 



ical Fermi surface and used an electron effective mass 
m* = 1.24mo, where mo is the bare electron mass, ex- 
tracted from Shubnikov-de Haas oscillations .7 1. In the 
d = 5.5 nm sample, we found £ rj 100 nm, £bcs ~ 470 
nm, Tt r ~ 6.2 x 10~ 14 s. Therefore, we conclude that our 
system is in the dirty limit, and spin-orbit coupling can 
be treated as a perturbation. This however becomes less 
clear with decreasing thickness, for which T tr /T so « 0.1. 
As a further check, we used a non-perturbative theory 
proposed by Schopohl et al. [H] 



The value of Hj! 2 for 



various A so using this theory is shown in Fig. S2 along 
with the original WHH model (neglecting the finite size 
corrections of Aoi et al. [5]). Since in absolute terms the 
observed values of i?| 2 are not large due to the relatively 
low T c , we estimate an error of only ~ 20 % in the de- 
termination of r so by using the original WHH model, as 
shown, which does not significantly affect the result of 
the original fitting. 
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